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ps| , Abstract 

^ , The intrinsic BL Lac spectra above few hundreds GeV can be very different from the 

I observed ones due to the absorption effects by the diffuse Extragalactic Background 

Light (EBL), at present poorly known. With the recent results, there are now 4 
sources with good spectral information: Mkn 421 (z=0.031), Mkn 501 (z=0.034), 
lES 1426+428 (z=0.129) and lES 1959+650 (z=0.047). Making simple assumptions 
on the shape of the intrinsic spectra (according to the present blazar knowledge), 
CN . we have considered the effects of different EBL spectral energy distributions (SED) 

\ for the first time on all 4 objects together, deriving constraints for the EBL fluxes. 

■ These resulted significantly lower than many direct estimates. 
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1 Introduction 

The interpretation of blazars TeV spectra is at present affected by a fun- 
damental ambiguity, due to the absorption of high energy photons by the 
'So . diffuse extragalactic background light (EBL) through 7-7 collision and pair 

production. The still large uncertainties on EBL knowledge prevent a reliable 
reconstruction of the intrinsic source spectrum, while conversely data from the 
^ ' synchrotron peak alone are not sufficient to univocally constrain the high en- 

■ ergy peak of the blazars Spectral Energy Distribution (SED), thus hindering 

the possibility to disentangle absorption from intrinsic features. Until more 
and better data are available, two possible lines of study are: 1) to assume 
a "best estimate" for the EBL and then to study the reconstructed intrinsic 
spectrum in the context of blazars emission models (see e.g. (j^); or 2) to make 
simple, "reasonable" assumptions on the shapes of the intrinsic TeV spectra 
(according to our present blazar knowledge), and then to limit the range of 
possible EBL fluxes (and spectral shapes) compatible with them. 

This second approach becomes more constraining as more sources are avail- 
able, since the EBL must be the same for all and how much an incident spec- 
trum is deformed depends strongly on redshift. With the recent results on lES 
1426+428 (3 IS E3) and lES 1959+650 (0; IB), we have now 4 objects with 
measured TeV spectra (see Fig. 1) and at 3 very different redshifts (~0.03, 
0.047, 0.129), thus allowing to draw multiple constraints. In the following, we 
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Fig. 1. Lower panels: EBL data points and the overall shape (continuous line) used 
to calculate the absorption-corrected TeV spectra showed in the corresponding up- 
per panels. Upper panels: filled symbols, the observed data; open symbols, the 
absorption corrected values. Circles: HEGRA data; triangles: CAT data; squares: 
WHIPPLE data. For lES 1959-^650, the CAT results are presented as the shape 
of the fitting model (a log-parabolic curve; details in (II)); the CAT and HEGRA 
data correspond to the high states in May-June and May-|-July 2002, respectively. 
As reference for the eye, in the lES 1959-1-650 and lES 1426-1-428 subpanels, two 
power-law spectra are also shown, with energy index a = 0.0 and 0.5, and 0.5 
and 0.9 respectively {Fi, oc z/~"). Left: EBL shape matching the reported measure- 
ments. Right: as left, but with lower fluxes in the 30-100 ^m range to avoid Mkn 
501 and Mkn 421 high energy upturns (as a reminder, the threshold condition for 
7-7 absorption is given by A < Xthr — 4.75(i?/lTeV) nm ). 

will focus on this second approach, and present the first preliminary results 
on the EBL SED obtained considering the effects on all 4 TeV sources at the 
same time. 



2 Data and Results 

Fig. 1 and 2 show different EBL shapes (obtained simply by connecting flux 
values at fixed frequencies) and the corresponding effects on the absorption- 
corrected TeV spectra. The compilation of measurements and lower limits 
have been taken from (Q), updated with the most recent publications fref. 
in (j^). For the blazar spectra, we have used the HEGRA data 111; Q; S S) 
conmleted at low energy by the CAT and WHIPPLE results for lES 1426-^428 
1^ [l3 ) and lES 1959-1-650 (11). In calculating the optical depths, we have 
assumed no model for EBL evolution, but only the changes in energy and 
photon density due to cosmology. For the latter, we adopted the recent W- 
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Fig. 2. Left: as Fig. 1 right, but with lower fluxes at 1-3 /im to avoid strange spectral 
shapes and spectra harder than a = 0. Right: further changing the overall shape 
to obtain an intrinsic lES 1959+650 spectrum between 0.5 and 1 TeV not harder 
than a ~ 0.5, and the lES 1426+428 spectrum similar to the BeppoSAX X-ray one 
(a ~ 0.9). 

MAP results, namely: Hq = 70 km/s/Mpc, Qmat = 0.3 and = 0.7. 



Our main assumptions for a "reasonable" blazar TeV spectrum have been: a 
smooth (power-law or curved) continuum, with no sharp features as up-turns 
or strong breaks, and a spectral index a {F^ cx z/~") not too hard, i.e. possibly 
> 0.5, which is the hardest value expected from an electron distribution cooled 
by synchrotron or inverse Compton mechanisms. Fig. 1 (left) show the results 
with an EBL shape designed to match the reported EBL measurements, in 
particular the high points at 1-3 yum (0t), and the first estimates at 60 /im 
(0). This gives origin to a sharp up-turn in the Mkn 501 and Mkn 421 spectra, 
and to an extremely hard spectrum up to ~1 TeV for lES 1959+650 and 
lES 1426+428 (a ~ —1, i.e. vF,j oc z/^), with a sharp break for the latter. To 
avoid these effects, it is necessary to lower the EBL flux in the range 30-80 /im 
(to avoid the up-turns in the two Markarians, Fig. 1 right) and in the range 
1-3 yum (to reduce the hardening of the low energy spectra, Fig. 2 left). Note 
also that the slope between 2 and 10 /xm should be around (or steeper than) 
\Fx oc A~^, otherwise the flattening of the observed lES 1426+428 spectrum 
above 2 TeV would imply again a sharp up-turn and very hard slope for the 
intrinsic one |l(] : 0) . 

With an EBL SED given in Fig. 2 (left), which is similar to some model 
predictions |lH), one can avoid strange shapes for the blazar intrinsic spectra, 
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but still has to face two uncommon features: a very hard spectrum for lES 
1959+650 {a ^ between 0.5 and ~2 TeV) and the overall SED of lES 
1426+428, whose TeV spectrum would imply a high energy peak at or above 
~8 TeV with a Compton dominance Lc/Ls > 10, unexpected in this class of 
objects (j^. To avoid this, the EBL fluxes should be of even lower intensity 
below ~3 /im , approaching the lower limits given by the HST counts 
With a shape as given in Fig. 2 (right), the intrinsic lES 1959+650 spectrum 
would have a low energy slope around a ~0.5, and lES 1426+428 would be 
characterized by a lower (more usual) Compton dominance (Lc^Ls) with a 
TeV spectral slope similar to the X-ray one (a;~0.9, (|^ Eoh). These limits can 
be loosened a bit taking into account the residual uncertainty on the absolute 
flux and energy calibration of Cerenkov telescopes, estimated in ~20% and 
~15% respectively (i.e. systematically decreasing energy and normalization of 
the data points of all the 4 TeV sources). 

3 Conclusions 

To avoid "unreasonable" shapes for the blazar intrinsic spectra, at least ac- 
cording to our present knowledge of the blazar phenomenon, the data of the 4 
TeV BL Lacs considered together put some limits to the EBL fluxes and spec- 
tral shape, which can be roughly given as: <15 nW/m^/sr in the range 0.8-3 
/im , <6 nW/m^/sr at 60yum (these values become <20 and <10 nW/m^/sr, re- 
spectively, including the systematic error), and with a slope around or steeper 
than XFx oc between 2 and 10 /xm . Conversely, if higher values or dif- 
ferent shapes should be confirmed by more constraining data or models, this 
would require a deep revision of our understanding of the blazars emission 
mechanisms and physical conditions. 
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